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ABSTRACT 

We construct both aligned and unaligned (logarithmic spiral) stationary configurations 
of nonaxisymmetric magnetohydrodynamic (MHD) disks fi'om either a full or a partial 
razor-thin power-law axisymmetric magnetized singular isothermal disk (MSID) that 
is embedded with a coplanar azimuthal magnetic field Bg of a non-force-free radial 
scaling r~^^^ and that rotates differentially with a fiat rotation curve of speed aD, 
where a is the isothermal sound speed and D is the dimensionless rotation parameter. 
Analytical solutions and stability criteria for determining are derived. For aligned 
nonaxisymmetric MSIDs, eccentric to = 1 displacements may occur at arbitrary 
in a full MSID but are allowed only with a^D^ ^ (7^/2 in a partial MSID {Ca is the 
Alfven speed), while each case of |to| > 1 gives two possible values of for purely 
azimuthal propagations of fast and slow MHD density waves (FMDWs and SMDWs) 
that appear stationary in an inertial frame of reference. For disk galaxies modeled by 
a partial MSID resulting from a massive dark-matter halo with a fiat rotation curve 



and a^D^ > 



stationary aligned perturbations of to = 1 are not allowed. For 



unaligned logarithmic spiral MSIDs with |m| > 1, there exist again two values of D^, 
corresponding to FMDWs and SMDWs that propagate in both radial and azimuthal 
directions relative to the MSID and that appear stationary in an inertial frame of 
reference. The larger is always physically valid, while the smaller is valid only 
for a > Ca/'2 with a positive surface mass density Sq- For observational diagnostics, 
we examine the spatial phase relationships among enhancements of gas density and 
magnetic field as well as velocity perturbations. These results are useful for probing 
magnetized bars, or lopsided, normal, and barred spiral galaxies as well as for testing 
numerical MHD codes. In the case of NGC 6946, interlaced optical and magnetic field 
spiral patterns of SMDWs can persist in a disk of flat rotation curve. Theoretical issues 
regarding the modal formalism and the MSID perspective are also discussed. 

Key words: accretion, accretion disks — galaxies: evolution — galaxies: magnetic 
fields — galaxies: barred, spiral — ISM: magnetic fields — MHD 



1 INTRODUCTION 

It is a challenge to study the large-scale dynamics of vari- 
ous morphologies of disk galaxies such as bars and lopsided, 
barred, and normal spiral structures (e.g., Baldwin, Lynden- 
Bell, & Sancisi 1980; Richter & Sancisi 1994; Rix & Zarit- 
sky 1995). To include gravitational interactions and magne- 
tohydrodynamics (MHD) of the interstellar medium (ISM) 
and magnetic field, we risk in making the task even more 
formidable. However, multi-wavelength observational diag- 
nostics involving both ISM and magnetic field do provide 



indispensable clues to the overall dynamics. For example, 
there have been growing numbers of high-quality observa- 
tions on large-scale magnetic field structures in nearby spi- 
ral galaxies (e.g., Sofue et al. 1986; Kronberg 1994; Beck et 
al. 1996; Beck 2001 and references therein). We here venture 
to formulate a limited yet nontrivial theoretical MHD disk 
problem in which stationary nonaxisymmetric MHD per- 
turbation configurations are constructed from a background 
axisymmetric MSID of interstellar gas medium embedded 
with a coplanar azimuthal magnetic field. We search for 
both aligned and unaligned stationary configurations in a 
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two-dimensional self-gravitating MSID with a flat rotation 
curve. 

Our model analysis here will give various stationary 
morphologies in a magnetized gas disk, including bars and 
lopsided, barred, and normal spiral structures. Moreover, we 
provide phase relationships of spatial patterns among mag- 
netic field, gas density, and velocity perturbations that can 
be examined observationally. Regarding the recent wavelet 
analysis on multi-wavelength data of the spiral galaxy NGC 
6946 (Frick et al. 2000, 2001) that revealed an extension of 
the interlaced magnetic and optical spiral structures into the 
outer disk with a largely flat rotation curve, our analysis here 
conveys an important message that stationary logarithmic 
spiral patterns of slow MHD density waves (SMDWs; Fan & 
Lou 1996; Lou & Fan 1998a), with interlaced spiral enhance- 
ments of magnetic field and gas density, can indeed persist 
in an extended MSID with a flat rotation curve (Lou & Fan 
2002). 

Conceptually, the stationary pattern problem here is 
closely tied to the density wave problem in a differentially 
rotating disk (Syer & Tremaine 1996; Shu et al. 2000). It 
is the balance between relevant wave pattern speeds and 
disk rotation that leads to possible nonaxisymmetric pat- 
terns stationary in an inertial frame of reference. Histor- 
ically, the seminal idea of MHD density waves was con- 
temporary (Lynden-Bell 1966; Roberts & Yuan 1970) with 
the early development of density wave theory four decades 
ago (Lin & Shu 1964, 1966; Goldreich & Lynden-Bell 1965; 
Toomre 1969, 1977; Lin 1967, 1987; Shu 1970a, b). While 
the stellar disk provides a massive "template" , MHD density 
wave processes in the ISM disk (Fan & Lou 1996, 1997, 1999; 
Lou & Fan 1997, 1998a, b, 2000a, b, 2002; Lou, Han, & Fan 
1999; LYF 2001; LYFL 2001) do have additional dynamic 
freedoms. 

Theoretically, there exist complemetary perspectives, 
different motivations, and independent approaches to study 
various bar phenomena in rotating self-gravitating fluid bod- 
ies. The masterpiece of Chandrasekhar (1969) on ellipsoidal 
figures of equilibrium summarizes the beautiful mathemat- 
ical descriptions of Maclaurin, Jacobi, Dedekind, Riemann 
ellipsoids and Poincare pear-shaped configurations as well 
as their close interrelations on the basis of instability analy- 
ses and characteristics of bifurcations in incompressible self- 
gravitating fluids. Complementary numerical studies of com- 
pressible, self-gravitating, differentially rotating, nonspher- 
ical equilibrium figures have shown striking family resem- 
blances to these classical exactly solved problems (Ostriker 
1978). By collapsing the dimension along the rotation axis, 
one can study equilibria and stability properties of two- 
dimensional Riemann disks of uniform rotation relevant to 
central parts of disk galaxies (Weinberg & Tremaine 1983; 
Weinberg 1983). Another broad class of disk problems in- 
volves stability properties of singular isothermal disks (SIDs) 
(Zang 1976; Toomre 1977; Lemos et al. 1991; Lynden-BeU 
& Lemos 1993; Goodman & Evans 1999). Syer & Tremaine 
(1996) found solutions to a class of stationary nonaxisym- 
metric perturbation SID configurations. Shu et al. (2000) 



derived solutions for stationary perturbation configurations 
in isopedically magnetized SIDs and interpreted them as on- 
sets of bar-type and barred-spiral instabilities (Galli et al. 
2001). Our analysis here parallels that of Shu et al. (2000) 
but with a coplanar magnetic field in an MSID, and we ex- 
amine the MSID problem from the perspective of stationary 
fast and slow MHD density waves (FMDWs and SMDWs; 
Fan & Lou 1996; Lou & Fan 1998a; LYF 2001). We derive 
a form of magnetic virial theorem for an MSID and suggest 
the ratio of rotation energy to the sum of gravitational and 
magnetic energies to be crucial for the MSID stability. 

Shu et al. (2000) pursued the analogies of bar-type in- 
stabilities known in thin self-gravitating disks (Hohl 1971; 
Miller et al. 1970) as well as in rotating self-gravitating 
spheroids and ellipsoids of uniform density (Chandrasekhar 
1969; Ostriker 1978) for equilibria and instabilities of isope- 
dically magnetized SIDs. In spite of the known difficulties 
(Galli et al. 2001), the fascinating scenario of the so-called 
"fission theory" by Poincare (1885), Liapunov (1905), and 
Jeans (1928) that regards incompressible ellipsoidal figures 
of equilibrium as potential candidates for fissioning through 
pear-shaped equilibria into binary stars has lead Shu et 
al. (2000) to speculate that nonaxisymmetric, magnetized, 
compressible, and perhaps truncated SIDs might prove to 
be promising candidates for fragmentation into binary- and 
multiple-star systems. They also hypothesized that a rapid 
loss of magnetic flux at a certain stage might hold the crucial 
key of resolving or bypassing the relevant known difficulties. 
In line with these speculations, the present MSID problem 
should be relevant to star formation research as well. 

Based on theoretical analyses (Lau & Berlin 1978; Lin & 
Lau 1979; Berlin & Mark 1979) to retain higher-order effects 
for the WKBJ expansion of the Poisson equation, Berlin et 
al. (1989a, b) heuristically argued and invoked the cubic dis- 
persion relation of density waves in a thin fluid stellar disk as 
the conceptual basis for classifying barred and normal spiral 
galaxies (viz., the Hubble classification scheme of galaxies; 
e.g., Lin & Lau 1979) constructed numerically by solving 
the standard linear integro-differential density wave equa- 
tions with chosen boundary conditions. Important aspects 
of this modal perspective have been comprehensively sum- 
marized in Berlin & Lin (1996). In their scenario, bars and 
barred spirals are essentially viewed as density-wave phe- 
nomena that are better described when the effects of long- 
range self-gravity and differential rotation are more fully in- 
cluded. Specifically, the cubic dispersion relation becomes 
cubic in the radial wavenumber k by including the tangential 
shear force (TSF) for nonaxisymmetric coplanar perturba- 
tions (LYF 2001). Besides the familiar short- and long- wave 
branches (both are somewhat modified by the TSF) , there is 
now a third wave branch of "open modes" characterized by 
even smaller k that bear striking resemblances to bars and 
barred spirals when superposed with an axisymmetric bulge 
(see Fig. 14 of Berlin et al. 1989a; Lin 1996 private commu- 
nications). As we will solve the Poisson integral exactly, it 
would be natural to pursue the correspondence between the 
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modal perspective of Bertin & Lin (1996) and the (M)SID 
results here as well as of Shu et al. 

This paper is structured as follows. We formulate in §2 
the problems of full and partial MSIDs and display relevant 
MHD equations. Solutions and analyses for aligned and un- 
aligned MSID configurations are presented in §3. In §4, we 
examine the spatial phase relationships between enhance- 
ments of gas density, magnetic field, and velocity distur- 
bances. In §5, we discuss the connection between the modal 
formalism and the MSID perspective for classifying normal 
and barred spiral galaxies, indicate the implication to multi- 
band observations of spiral galaxy NGC 6946, and summa- 
rize the results. Mathematical formulae are collected in Ap- 
pendices A— E for the convenience of reference. 



2 FORMULATION OF THE MSID PROBLEM 

The key difference of our formulation and that of Shu et 
al. (2000) is that their magnetic field is poloidal thread- 
ing across the disk and may be effectively relegated into 
two dimensionless parameters* © and e (Schmitz 1987; Shu 
& Li 1997) such that Qa^ and G ~* eG {a is the 

isothermal sound speed and G is the gravitational constant), 
whereas the non-force- free magnetic field of our model is az- 
imuthal and coplanar with the disk. In essence, the effect of 
their poloidal magnetic field can be scaled away such that 
the analysis is equivalent to a hydrodynamic one, while the 
coplanar magnetic field of our model can give rise to new 
classes or features of stationary nonaxisymmetric MSID con- 
figurations. We refer to our model as an MSID problem, al- 
though the SID problem as formulated by Shu et al. (2000) 
is isopedically magnetized, and proceed to examine several 
specific aspects in an orderly manner. 



2.1 Nonlinear unsteady MHD equations 

For large-scale MHD processes, dissipative effects may be 
ignored as a first approximation. In cylindrical coordinates 
{r,0, z), ideal time-dependent MHD equations include 

1 dm 



as 1 djrEu) 
dt r dr 



+ 



= 



2.1.1 



for the mass conservation, where E is the vertically inte- 
grated mass density, it is the radial velocity, j = rv is the 
specific angular momentum in the vertical z direction and v 
is the azimuthal velocity in ^—direction. 



du 



du 



j du f 



1 9n 



dt 



dr 



dr 



dzBo 
Anr 



E dr 
dirBg) 
dr 



dBr 

de 



(2.1.2) 



* Specifically, 6 = (A^ -f 1 + 2»y2)/(A2 + rf) and e = 1 - A^^ 
with A = 2TrG^/^T,/B^ = const, and 'q = Igy |/(27rGi;) where 
is the coplanar gravitational acceleration. 



for the radial momentum equation, where H is the vertically 
integrated gas pressure and (j>T is the total gravitational po- 
tential inclusive of that from a possible dark-matter halo. 



I i dj ^ ^ 9n 

dt dr r^de T, de 



d(j>T 1 

~de' s 



dzB-r 
4-K 



djrBe 
dr 



dBr 



(2.1.3) 



for the azimuthal momentum equation in ^—direction 

s(C,^,t)CrfC 



= -G ^ dil^J^ 



(2.1.4) 



[C2 + r2-2Crcos(V'-6i)]i/2 

for the Poisson integral equation where F<I)t is the gravita- 
tional potential from the gas disk with < F < 1 , 
d{rBr) dBg 



dr 



= 



(2.1.5) 



for the divergence- free condition on _B = [Br, Bg,0), 

-rn- = ZT^i^Bg^vBr) 



(2.1.6) 



-{uBg — vBr) 



(2.1.7) 



for the radial magnetic field induction equation, and 
dBe _ d_ 

dt ^ dr' 

for the azimuthal magnetic field induction equation. The 
polytropic approximation H = a^E is invoked as usual where 
a is the isothermal sound speed. 

For a static vertical balance without a vertical fiow ve- 
locity Vz, the vertical force balance takes the form of 







dz / dz 



dzli^l±^ , (2.1.; 



dz Stt 

where gas and magnetic pressures work together against the 
vertical gravity towards the disk ai z — 0. 

2.2 Equilibrium of an axisymmetric MSID 

For an MSID geometry of axisymmetry, one solves the MHD 
equations with the Poisson integral (2.1.4). We presume a 
fiat rotation curve with a background disk angular rotation 
speed fl — aD/r , where D is the dimensionless rotation 
parameter. The cases of D > 1 and D < 1 correspond then 
to supersonic and subsonic MSID rotations. To avoid the 
magnetic field winding dilemma (Lou & Fan 1998a), the 
background magnetic field is taken to be azimuthal about 
the symmetry axis with a scaling of Be = J-r~^^^ where 
is a constant. The Lorentz force of this Bg profile ^ is 
included in the MSID radial force balance. The vertically 
integrated gas pressure Ho of the background is related to 
the background surface mass density Eo by Ho = a^Eo. The 
epicyclic oscillation frequency k of a flat rotation curve is 
defined by = {2n/r)d{r'^Q,)/dr = 2n^ and the Alfven 
speed Ca in an MSID is defined by Cl = / d2B|/(47rEo). 

The rotation speed of a disk galaxy Vg = Q,r — aD ^ 
Ca due to the presence of a massive dark-matter halo (pre- 
sumed to be axisymmetric here). By the equipartition argu- 
ment, the thermal and magnetic energy densities are compa- 

t Bg scales as r~^ for a force-free azimuthal magnetic field. 
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Table with a ~ Ca - It then follows that D ^ 1 for supersonic 
and super-Alfvenic rotations in typical disk galaxies. 

To attribute a fraction (1 — _F) of the total gravity in 
the background equilibrium to an axisymmetric dark-matter 
halo that is unresponsive to gas disk perturbations, one may 
write dcpr/dr — Fdipr/dr + (1 — F)d4>T/dr . The case of 
_F = 1 is referred to as a full (M)SID, while the case of < 
F < 1 is referred to as a partial (M)SID (Syer & Tremaine 
1996; Shu et al. 2000). In reference to equations (5), (6), 
(21) and (27) of Syer & Tremaine (1996), their parameter 
/ is related to F here by / = (1 — F)/F. The background 
rotational equilibrium of an MSID requires 

g^D^ d4>T C\ 

r r dr 2r ' 

where the portion of the gravitational potential associated 
with the gas disk satisfies Fd(j>T/dr = 27rGEo in the disk 
plane at 2: = 0, as required by the Poisson equation in a 
razor-thin MSID. In force balance (2.2.1), the gravity and 
net Lorentz forces are radially inward (i.e., an outward mag- 
netic pressure force but a stronger inward magnetic tension 
force), while the gas pressure and centrifugal forces are ra- 
dially outward. It follows for an MSID that 



(2.2.1) 



So = F[a^{l + D^) - Ci/2]/(27rGr) 



(2.2.2) 



In the context of a protostellar disk, one cannot invoke a 
dark- matter halo and one thus deals with a full (M)SID of 
F = 1. Shu et al. (2000; also Galh et al. 2001) recently 
performed extensive analysis on constructions of stationary 
nonaxisymmetric perturbation solutions from a background 
axisymmetric SID that is isopedically magnetized with a 
poloidal magnetic field. These solutions can be classified as 
aligned and unaligned configurations (Kalnajs 1973) that re- 
quire specific values of for their very existence. These re- 
sults in idealized theoretical settings are interesting and im- 
portant in that these solutions, with proper interpretations, 
may be pertinent to configurations of bars, barred spirals, 
and normal spirals in a differentially rotating disk, and may 
bear implications to the onset of bar or bar-type instabili- 
ties that lead to formation of these possible configurations. 
For example, for the aligned case of nonaxisymmetric sta- 
tionary perturbations (of angular variation exp[— imS] with 
an integer m) studied by Shu et al. (2000) (see their eqns 
[25] — [27]), one finds an unconstrained or arbitrary for 
|m| = 1 (see also equation [46] of Syer & Tremaine 1996), 
andD^ = |m|/(|m| -t-2) for |m| > 2. The latter with < 1 
implies a subsonic protostellar disk rotation. 



2.3 Coplanar MHD perturbations in an MSID 

From equations (2.1.1) — (2.1.7), one derives coplanar MHD 
perturbation equations in an MSID of gas. 



dui dui 2nji 



d ( a^Ei 



dr 



" 2Eor 



_1_ r dzBg djrbg 
Eo J inr dr 

_1_ f dzbe djrBg) 
Eo / 47rr dr 



dbr 



(2.3.2) 



dt ^ 20. ^^ de ^ de 



Eo 

dzbr dirBg) 
4:TT dr ' 

Ei(C,^,t)CdC 



[C2-Fr2-2Crcos(V'-6')]i/2 



dbr 

'dt 



djrbr) 
dr 

- 

rde 
d 



dbg_ 

de 



= 



(2.3.3) 
(2.3.4) 
(2.3.5) 



(uiBg — rflbr) 
(uiBg — rQ,br) 



(2.3.6) 
(2.3.7) 



dbg_ _ 
dt dr ^ 

where b = {br,bg,0) is the coplanar magnetic field pertur- 
bation, is the perturbation of F<j)T from gas distribution, 
and other variables with subscript 1 are perturbations to the 
pertinent equilibrium variables. Here, we set — b^ = and 
do not consider vertical variations across the disk. All with 
the harmonic exp{iujt — im9) dependence, we introduce com- 
plex radial variations S{r), U{r), J{r), V{r), R{r), and Z{r) 
for El, Ml, ji, 01, br, and bz. Coplanar MHD perturbation 
equations (2.3.1) — (2.3.7) can then be reduced to 

zmEo 



i(Lj - mn)S + --^(rSoU) ■ 
r dr 



-J = 



(2.3.8) 



i{uj — mQ)U - 

C\S 
+-. 



20/ 

r 

1 

" e;^ 



9$ 

dr 
dzBg 



2Eor Eo / 47rr 
where $ = a^S/Eo + V 



_i_ f dzZ djrBg) 
Eo J A-Kr dr 

d{rZ) 



dr 



i{Lj — m.Q,)J + -^^U 



1 



imR 



dzR d(rBg) 



V{r) ^-G (pdx 



Eo / 47r dr 

S{C)cos{mx)CdC 
g +r'^ -2CrcosxV^^ 

d{rR) 



dr 



imZ = 



./ T., imBg 
i[uj - mQ.)RA -U ■ 



iLoZ = ^(rVLR) - -^(BgU) 



(2.3.9) 

(2.3.10) 

(2.3.11) 

(2.3.12) 
(2.3.13) 
(2.3.14) 



dr dr 

It suffices to use only two of the three equations (2.3.12) — 
(2.3.14) for magnetic field perturbation b = (6,.,?>e,0). 

Rearrangement of equations (2.3.8) — (2.3.14) in terms 
of S, V, J and iU further leads to 



Z ^- 



i d{rR) 
m dr 



(2.3.15) 
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from the divergence-free condition (2.3.12) of b, and 

R = (2.3.16) 
r(u — mil) 

from the radial magnetic induction equation (2.3.13). Using 
equations (2.3.15) and (2.3.16) in the radial and azimuthal 
momentum equations (2.3.9) and (2.3.10), one derives 

i{uj - mQ.)U = + 

r or 2Eor 



1 d 



d 



^1/2 Qj. I Qj. 

1 d 



u 



+ 



r-^(Lj — mfi) 2r^/'^ dr 



rV2(a;_mn) 
U 



and 



i{uj - mfl)J + —U = im$ - — ^— — 

2S2 2r(a; — mil) 



(2.3.17) 



(2.3.18) 



2.4 Stationary MHD perturbations in an MSID 

For nonaxisymmetric MHD perturbations stationary in an 
inertial frame of reference, we set a; = in equations (2.3.8), 
(2.3.16) - (2.3.18) and obtain 

i^(.E„z[/) + !^J = 0, (2.4.1) 
r or 



mns 

2Q.J 



mQiU H = -TT- 



r 

milj 4 



dr 



CIS 
2Eor 

iU 

dr \ mO.r^/'^ 



d 



iU 



d 



2^1/2 Qj. \ ^{7^1/2 

+ -^T^ , (2.4.2) 



rK.'^ ^ CliU 



iR ■- 



BgiU 

fir 



(2.4.3) 



(2.4.4) 



together with equations (2.3.11), (2.3.15), and the definition 
of $ after equation (2.3.9). As the pattern speed uip = uj /m 
of possible MHD density waves is set to zero in an iner- 
tial frame of reference a priori, these equations are to be 
solved to determine the proper values of the dimensionless 
rotation parameter D. As the system may support possible 
FMDWs and SMDWs (Fan & Lou 1996; Lou & Fan 1998a), 
two proper values of D are expected. 



3 SOLUTION ANALYSIS OF STATIONARY 
MSID CONFIGURATIONS 

3.1 The case of aligned perturbations 

The case of jn = should be handled with care. It would be 
misleading to use equations (2.4.3) and (2.4.4). One should 
examine equations (2.3.8) — (2.3.14) with ui = m = 0. By 
setting U = R = 0, equations (2.3.8), (2.3.10), (2.3.12) 
and (2.3.14) are satisfied, and equation (2.3.13) is identi- 
cally zero. There is no constraint on Z. However, by setting 
Z (X r~^^'^, S ix r~^, Jcx r, V oc Inr, and $ oc const. -I- In r. 



remaining equations (2.3.9) and (2.3.11) can be made consis- 
tent with a rescaling of the axisymmetric background MSID. 

For the nonaxisymmetric aligned case with m 7^ 0, one 
takes the density-potential pair 

rS' = const. (3.1.1) 

and 

1/ = -(27rG/|?n|)r5' = const. (3.1.2) 

(e.g., Shu et al. 2000; Galh et al. 2001). For a constant iU 
(see eq. [3.1.4] below), the mass conservation (2.4.1) requires 



(3.1.3) 



and the azimuthal momentum equation (2.4.3) gives consis- 
tently a constant ill by 

iU = m($ - nV^S'/Eo)/[Ci/(2J7r) - Qr] . (3.1.4) 

By equations (3.1.1), (3.1.2), and (2.2.2), it follows that $ = 
a'^ S/T,o + V = const, and a combination of equations (3.1.3) 
and (2.4.2) leads to 



mQr 



Cl{m^ - 1/2) 
mQ,r 



Substitutions of expressions (3.1.2) — (3.1.4) into equation 
(3.1.5) give the solution condition or stationary dispersion 
relation for aligned nonaxisymmetric MHD perturbations, 

m^nr - Clim'^ - l/2)/{nr) f _ 2tvG _ »V 
Cl/{2nr) - fir l^E^ ^ ~\m\ e7 

Eo 2Eor 

For a later examination of the spatial phase relationship 
between bg and Ei, we further derive the following results. 
By equations (2.3.15) and (2.4.4), the 0- component of the 
magnetic field perturbation is related to ill by 

Z = -iUBe/{2mnr) . (3.1.7) 

By using expression (3.1.4) for constant ill and expression 
of constant <1> in equation (3.1.7), one obtains 

BeS (a^ - 27rGEor/|m| - nV^) 



2Eo {Cl/2 - 02r2) 

The solution criterion (3.1.6) may then be written as 

2nV2 - Ca/2 



(3.1.8) 



27rG'Eor/jmj - fiV^ 
Cl/2 - n^ri 



m^n^r^ - Cl(m'^ - 1/2) ' 
(3.1.9) 

By equation (3.1.8), the sign of either left-hand side (LHS) 
or right-hand side (RHS) of equation (3.1.9) will determine 
the phase relation between the surface mass density pertur- 
bation S and the azimuthal magnetic field perturbation Z. 

3.2 The physical nature of the solution criterion 

To fully understand the nature of solution condition (3.1.6), 
we examine two conceptually related cases in order. The first 
case of Ca = is essentially the one studied by Shu et al. 
(2000) (e = 1 and O = 1) and can be explained as purely 
azimuthal propagation of hydrodynamic density waves (Lin 
& Shu 1964, 1966). By this clue, we proceed to show that the 
second case of 7^ corresponds to two possible situations 
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of purely azimuthal propagations of FMDWs and SMDWs 
(Lou & Fan 1998a). 

3.2.1 The aligned full and partial SID cases 

We have verified with = and _F" = 1 in our analy- 
sis that the resuhs (25) - (27) of Shu et al. (2000) come 
out naturally, namely, either is arbitrary for \m\ = 1 or 

= \m\/{\m\ + 2) for |m| > 2. The latter with < 1 
corresponds to a subsonic SID rotation. 

Moreover, the solution condition (3.1.6) or (3.1.9), with 
C\ — 0, can be written in the informative form of 



rn^r^^ = 20,'^ + rnajr'^ - 27rGEo|m|/r 



(3.2.1) 



In reference to the well-known WKBJ dispersion relation 
of density waves (Lin & Shu 1964, 1966 or equation [39] of 
Shu et al. 2000), equation (3.2.1) can be readily obtained 
by replacing the radial wavenumber \k\ with the azimuthal 
wavenumber |m|/r and setting cij = in an inertial frame of 
reference. This clearly describes an azimuthal propagation of 
hydrodynamic density waves, and a stationary pattern in the 
sidereal frame of reference requires specific values of for 
different |m| values. It was pointed out by Shu et al. (2000) 
that perturbations of the aligned case has no radial wave 
propagation. We here provide a transparent physical inter- 
pretation for the aligned case of nonaxisymmetric stationary 
perturbations in terms of a purely azimuthal density wave 
propagation (retrograde relative to the disk) advected by the 
SID rotation such that the stationarity is sustained. As the 
corotation is at infinity for a stationary perturbation and 
— 257^, a solution of equation (3.2.1) appears outside the 
Lindblad resonances also located at infinity. Specifically, the 
solution is valid inside the inner Lindblad resonance (ILR) in 
a finite radial range. Formally, only the solution of |m| = 1 
appears within the Lindblad resonances. 

Equation (3.2.1) is quadratic in |m| and the two pos- 
sible values of |7ti| should be reminiscent of the long- and 
short-branches of density waves in terms of the radial wave- 
length 27r/|fc| (e.g., Binney & Tremaine 1987), even though 
the rn}V? term on the LHS of equation (3.2.1) is now in- 
volved in determining the proper values of |m|. By this per- 
spective and for a full SID with F = \, one might view the 
two values of jmj given by equation (3.2.1), namely |m| = 1 
with an arbitrary value, and |m| = 2D^/(1 — D^) with 
|m| > 2, as "long" and "short" stationary azimuthally prop- 
agating density waves. By the latter, one requires < 1 
and in the limit of — > 1, the value of |m| increases to 
infinity corresponding to extremely short azimuthal wave- 
lengths. Perhaps, the |m| = 2 case is the most interesting 
one that mimics a stationary barred configuration. 

For a partial SID with < F < 1 and C\ = 0, criterion 
(3.2.1) becomes 

rnD^ = 2D^ +ni - \m\Fijf -\- 1) . (3.2.2) 

The two values of |m| are then given by 

_ -F(p^ + 1) ± {F'^jlf- + \f -I- 8D^(I)^ 
" 2(D2 _ 1) 



11/2 



(3.2.3) 



where negative values of |m| must be rejected. 

To examine the solution property of (3.2.3) for |m|, con- 
dition (3.2.2) may be cast into the revealing form of 

D — |m| (|m \- F)l(jn \ F\m\-2) . (3.2.4) 

For a full SID with F = 1 in equation (3.2.4), one would 
have = \m\/{\m\ -I- 2) < 1 for a subsonic SID rotation. 
For F = 1 and jmj = 1 in equation (3.2.4), the value of 

is arbitrary for a nontrivial eccentric solution (Shu et al. 
2000), including the case of a cold rotating disk of = 
but Qr = aD 7^ (see discussions of Syer & Tremaine 1996 
about their eq. [43]) and the case of a non-rotating disk of 

= (see discussions of Syer & Tremaine 1996 about their 
eq. [45]). In the limit of F = 0, equation (3.2.4) gives 



= rnjirn - 2) > 1 



(3.2.5) 



(see eqns. [8], [21], [25], [44] of Syer & Tremaine 1996 with 
their /? = 0) for a supersonic SID rotation with jmj > 2, and 
the case of jmj = 1 is forbidden because > is required. 

Similarly for < F < 1, the case of jmj = 1 is no longer 
allowed for a partial SID because = — 1 is unphysical. In 
this context, we note that the corollary of Syer & Tremaine 
(1996) after their equation (44) is inaccurate in the sense 
that the situation of their /9 = (fiat rotation curve) and 
finite / 7^ (0 < F < 1) does not allow for aligned m = 
1 perturbation solutions. There are two distinct classes of 
solutions for jmj in general. For jmj < F~^, one would have 
a supersonic SID rotation with > 1, while for jmj > F~^, 
one would have a subsonic SID rotation with < 1. It is 
then clear that the case of F = 1 excludes the possibility of 
> 1, while the case of F = excludes the possibility of 
< 1. For < F < 1, aU solutions with jmj > 2 satisfying 
condition (3.2.2) appears inside the ILR. 

For spiral galaxies, a partial SID model with < F < 1 
is more relevant due to the presence of a massive dark-matter 
halo and there seems to be no need to worry about the 
jmj — 1 case of aligned eccentric displacements because they 
only exist for a full SID with F = 1 (Zang 1976; Toomre 
1977; Shu et al. 2000). However, in star-forming clouds not 
involving an axisymmetric dark-matter halo, such case of 
jmj = 1 for aligned eccentric displacements might lead to 
alternative SID configurations of interest (Shu et al. 2000). 



3.2.2 Aligned MSID perturbation configurations 

With C\ 7^ 0, the stationary criterion (3.1.6) or (3.1.9) for 
the aligned case may be cast into the informative form of 



m it 



2^f + 



C 



2 2 

A O, 

— H 



+- 



2-kGT.q 



27rGS o 

jmjr 

2 1 

' ~2 



4^2 



0, (3.2.6) 



which reminds us of the dispersion relation for spiral 
FMDWs and SMDWs (Fan & Lou 1996; Lou & Fan 1998a; 
LYF 2001) by simply replacing the radial wavenumber jfcj 
with the azimuthal wavenumber jmj/r and by setting lj = 
in an inertial frame of reference. In the earlier derivation 



© 2001 RAS, MNRAS OOO.ITII20I 



Stationary MSID configurations 7 



of WKBJ dispersion relation for FMDWs and SMDWs, we 
took a force-free Be oc . The present case of Bg cx r"^''^ 
is not force-free and some extra terms in expression (3.2.6) 
are due to the additional radial and azimuthal Lorentz force 
terms as well as geometric effects in equations (2.3.2) and 
(2.3.3). Physically, the two possible values contained in 
equation (3.2.6) should relate to the two situations of purely 
azimuthal propagation of FMDWs and SMDWs (retrograde 
relative to the disk) such that the disk rotation renders them 
stationary in an inertial frame of reference. 

The equivalent criteria (3.1.6), (3.1.9), and (3.2.6) may 
be arranged into the form of 



1/2)] 

X [a(l - D'^) - F{aD'^ +a - C\/2)/\m\] 



Cl/2) = 



(3.2.7) 



For a full MSID with F 
can be factored into 



1, the LHS of expression (3.2.7) 



{\m\-l){{a^D'' ~C\){l^D'')c 
+ (a2i32-Ci)(Ci/2- 
-[Cl/2-a\D^ + l)]Cl/2}^Q 



(3.2.8) 



Therefore, similar to the case of perturbations in an isope- 
dically magnetized full SID (Shu et al. 2000), the value 
is unconstrained for |m| = 1. For \m\ > 1, the second factor 
on the LHS of equation (3.2.8) may be rearranged into 



|m|(2 + lml)a^D^ 



[Ci/2 + 5CX\m\/2 + {a^ + Cl)\m\^]aD^ 
+ {Cl/2 - a^)Cl/2 + Ci\m\/2 + cW\m\^ = (3.2.9) 

with the last coefficient being positive for jmj > 1. By 
(3.2.9), there are two positive values of a?D^ , namely 

2 n2 _ ^^/2 + 5Ci|m|/2 + (a^ + Ci)|mp 



2|ml(2 - 



+ 



+ - 



2 2 

-4jm|(2 -I- |m|) 



+ {a^ + C\)\m\^ 



2 



1/2 



[2|m|(2 - 



(3.2.10) 



where the determinant A is non-negative for |m| > 1 (see 
Appendix A or equation [-D9]). 

For a partial MSID with < F < 1, expression (3.2.7) 
may be cast into the form of 

(2 - - F\m\)a^D'' + {A-B + C)aC)^ 

-'{C\lm')(m^ - l/2)A + Cl/4 = , (3.2.11) 

where A, B, C are defined by equations (Bl), {B2), {B3) 
in Appendix B. We emphasize that \m\ = 1 is no longer 
a solution with unconstrained for a partial MSID. One 
may solve equation (3.2.11) for two values of a^D^ 



aD' = {B-A-C± A'''^)/[2(2 - - F 



(3.2.12) 



where the determinant A is positive for \rn\ > 1 (see defini- 
tion [_B4] in Appendix B). 



For the case of |mj = 1, criterion (3.2.7) becomes 

(a^D^ - Ci/2)(1 - F)(a^D^ + - Cl/2) = , (3.2.13) 

where < F < 1. As 27rGrEo = F[a^{l + D'^)~Cl/2] > is 
a necessary requirement, the only possible solution of a^D^ 
for a partial MSID becomes a^D^ = (7^/2 . 

3.3 Onset of bar-type instabilities in MSIDs 

For the disk stability problem, there exists a key parameter 
for both secular and dynamic bar-type instabilities, namely, 
the ratio of the kinetic energy of rotation T to the abso- 
lute value of the gravitational potential energy W (Ostriker 
& Peebles 1973; Ostriker 1978). If one regards the station- 
ary configuration of aligned \m\ = 2 MSID as the onset of 
secular bar-type instabilities (Shu et al. 2000), it is then of 
interest to properly estimate the ratio T/|W|. For an MSID 
of infinite radial extent, both quantities T and W are infinite 
but their ratio remains finite. Starting from the radial force 
balance of the background magnetorotational equilibrium 

-Eo^V = -dn/dr - T.Qd<j>T/dr - T.oC'\/{2r) , (3.3.1) 

we multiply equation (3.3.1) by 2-Kr'^dr and integrate from 
to a radius R which is allowed to approach infinity eventu- 
ally. After an integration by parts for the gas pressure term 
up to R, one obtains the MSID virial theorem 

2{T + U)+W-M=2TTR^Ii{R) (3.3.2) 

within R, where by using equation (2.2.2) for Eo 



T; 



Eor^f2^27rr-dr ; 



a^D^F 
2G 



1 + D^ 



2a2 



R (3.3.3) 



is the kinetic energy of the MSID rotation. 



U = 



Il2nrdr — 



a'F 
G 



is the gas internal energy, 

yV=-[ r — Eo27rrdr = 
Jo 



a^F 
G 



1+D^ 



'2a2 



R (3.3.4) 



R (3.3.5) 



is the total gravitational potential energy (if there is no mat- 
ter outside R), 



M 



JdzB, 



8-K 



2G 



R (3.3.6) 



is the MSID magnetic energy within R. The usual key ratio 
is then 

T/\W\ = a^DV[2(a^ + aB)^ - C\l2)\ . (3.3.7) 

In view of the MHD virial theorem (3.3.2) in which W < 
by equation (3.3.5), one might suspect the quantity W — 
to play the role of W in the absence of coplanar magnetic 
field. Because 



a^(\ + D'^)F 



2 



R<0 (3.3.8) 



for Eo > 0, it would be suggestive of a new modified ratio 

T/\W-M\=D'^/[2{1 + D'^)] (3.3.9) 

to determine the instability criterion of an MSID for nonax- 
isymmetric aligned MHD perturbations. 
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3.4 The unaligned case of stationary 

nonaxisymmetric MSID configurations 

We now consider the unaligned case for logarithmic spiral 
structures (Kalnajs 1973; Shu et al. 2000). Combinations of 
azimuthal momentum equation (2.4.3) with radial momen- 
tum equation (2.4.2) and with mass conservation (2.4.1) give 



m it r — K r + Ca 



f c\ 



\ 2^27-2 



2n2r2 



■ raVlr — + 2mf2r$ — 

or 2Eo 



d f mO.rS m$ 
mQ.rS 



So 



m<I> 



(3.4.1) 



and 

mfirS d{iU) 



So 



dr 



2^72.2 2^.2, 'f/-S-0- (3-4.2) 



The Poisson integral (2.3.11) is a linear integral equation 
and has a number of complete sets of eigenfunctions, de- 
pending on the radial range of nonzero surface mass density 
(Snow 1952). For a radial domain of (0, 1), the radial part 
of the eigenfunctions may be Legendre polynomials (Hunter 
1963) or Bessel functions (Yabushita 1966), while for a ra- 
dial domain of (0, oo), one may use Bessel functions (Toomre 
1963) or logarithmic spirals (Kalnajs 1965, 1971). We here 
study unaligned logarithmic spiral structures. 

For logarithmic spirals (Kalnajs 1965, 1971) with con- 
stant pitch angles, the known potential-density pair is 

S = sr-^''2+'" and V = ^-1/^+'" , (3.4.3), (3.4.4) 

where a is a constant parameter that characterizes the radial 
variation of perturbations and the two constant coefficients 
s and V are related by 

V = ~2TiGNrr,{a)s , (3.4.5) 



with Nm{ct) = K{a,m) 

_ 1 r[(m + 1/2 -f ja)/2]r[(m + 1/2 - ia)/2\ 



(3.4.6) 



2 r[(m + 3/2 + ia)/2]V[{m + 3/2 - ia)/2\ 

being the Kalnajs function (Section IV of Kalnajs 1971) and 
r(...) being the gamma-function (Qian 1992; Appendix A of 
Syer & Tremaine 1996). The Kalnajs function is real and 
positive. Table 1 of Kalnajs (1971) contains numerical values 
of K[a,m) for different values of a and m. For large m or 
large a, one has Mm{,ci) = (m? + a^)~^^'^ approximately. For 
m = 0, 1, 2, the curves of K(a, m) versus a were displayed in 
Figure 1 of Shu et al. (2000). The Kalnajs function is even 
in m and a. With a sufficient accuracy (Shu et al. 2000), 
one may approximately take 

AA™(a) ?^ (m" +o? + 1/4)-'/" . (3.4.7) 

The exact recurrence or recursion relation for Nm{ct) is 

M^+i{a)M^{a) = [(m + 1/2)" + a"]"' (3.4.8) 



(Kalnajs 1971). Suppose one takes, for example, 

N2{c) ^ (4 + a"-hl/4)-'/" 

by approximation (3.4.7) with |m| = 2, it then follows suc- 
cessively from the recursion (3.4.8) that 

AAi(a)?^=(l-h 1/16 + aV4)'''V(l + 1/8 + aV2) < 1 
and 

^ (a2 + i/4)(l + 1/16 + a2/4)i/2 
(see Fig. 1 of Shu et al. 2000). To achieve a higher accuracy 
for Afm{a), one starts from approximation (3.4.7) with a 
large |mj and use the recursion (3.4.8) to successively derive 
Nm[oi) of descending |m|. 

Consistent with the logarithmic spiral forms of S{r) and 
Vij-) by equations (3.4.3) and (3.4.4), we take U to be 

iU = mr-'/"+'" , (3.4.9) 

where « is a constant coefficient. It follows from equations 
(2.4.1), (2.4.2), (3.4.5), (2.3.12), and (2.3.15) that 
mfls . f . , C\ 3^ ™ / ..2 



2Q?r'^ 2 



mas , 



2o2 2 2 2 



ci 



2Q?r'^ 



2n2r2 



= C\\ic + 2- 



+mQ.r{ia + ?,/2) 



rEo 



mQ,s 
s 

+ v 



.1/2 



iR ■- 



Be 



Z — — lur 

milr 



Qr V rEo 
2Eo 



BgiU 

laBe 



+ v 



iU 



(3.4.11) 



(3.4.12) 



(3.4.13) 



Equation (3.4.10) relates s and iu by relation (3.4.5). In- 
dependently, equation (3.4.11) is another relation for s and 
iu by relation (3.4.5). A combination of the two resulting 
relations in terms of s and iu then gives rise to the solu- 
tion criterion for stationary logarithmic spirals in an MSID. 
With an exact cancellation of the imaginary part, the solu- 
tion criterion or dispersion relation for stationary unaligned 
logarithmic spirals becomes 



rE( 



+ v 



milr 



2nr 
m 



2Eor V 2Q^r 



m 



rEr 



+ v 



a — 1 m 
milr ilr 



mils 



Eo 



(3.4.14) 



where s and v are related by equation (3.4.5). A physically 
more revealing form of criterion (3.4.14) is 

m'^Q.'^ - {20.'^ + \C\lr^ + a^ /r^ - 27rG'7V™(a)Eo/r] 
x{m^ + a +1/4) -2Ci/r"}m"f7" 
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+{m^C\/r^){[a'/r'^ - 27rGAA„(a)Eo/r] 

X (mVQ^ + 1/4- 1/2) --Ci/(4r^)} = . (3.4.15) 

In parallel, the two stationary solution criteria (3.2.6) and 
(3.4.15) for the aligned and unaligned cases correspond to 
each other remarkably well, especially in view of the effec- 
tive wavenumber {m? + + l/Af^r-^ (Shu et al. 2000). 
Again, the form of equation (3.4.15) reminds us of the dis- 
persion relation for spiral FMDWs and SMDWs (Lou 1996; 
Fan & Lou 1996; Lou & Fan 1998a; LYF 2001). These log- 
arithmic spiral MHD density waves propagate in both ra- 
dial and azimuthal directions relative to the MSID. For sta- 
tionary logarithmic spirals in an inertial frame of reference, 
criterion (3.4.15) leads to two possible values of a^D^. As 
both = and Eq contain the parameter, the de- 
termination of should be somewhat different from the 
standard WKBJ wave results derived in a background that 
is prescribed a priori. For example, for a full SID that is iso- 
pedically magnetized, Shu et al. (2000) obtains only one real 
solution of Of at a given for each m > 1 as shown in their 
Fig. 3 instead of the two branches of long- and short-waves 
(see also Syer & Tremaine 1996 for a discussion). While for 
a sufficiently large value in the special case of m = 0, 
there are two values of a as shown in their Fig. 2, bordering 
the ring fragmentation regime. 

Equation (3.4.15) may be cast into the form of 



3^ - - 27rGA/'m(Q)Eor] 
X [(m^ + + 1 /4)nV^ - (m^ + a 



1/4)C^ 



= -a^r^V* - (Ca/2 - 3nVV2)^ . (3.4.16) 

After a straightforward rearrangement using the Eq profile 
(2.2.2), criterion (3.4.16) becomes 

[m' + + 1/4 - (m^ + a' - l/4)Ci/(a'D^)] 
x{a' - [(1 + D'')a^ - Cl/2]FUm{a)} 

+ (Ci/2)[Ci/(2a"D") - 3/2] = . (3.4.17) 

For C\ = and F — 1, this criterion reduces to equation 
(37) of Shu et al. (2000). Also with C\ = Q and ^ in 
the so-called "breathing mode" regime, one obtains 

m^D^ = 2D^ + {m^ + l/4)[l~{l + D^)FM^{Q)] , (3.4.18) 
which differs from condition (3.2.2) derived earlier for the 
aligned case, especially in the regime of large |m|. 

In terms of a quadratic algebraic equation for a^D^, 
equation (3.4.17) may be written in the explicit form of 

[m^ - 2 + (m^ + l/4)FAA™(a)]a'*D'* (3.4.19) 



dim'' 



4 







where A, B, C are defined by equations (CI), (C2), (C3) in 
Appendix C. For |m| > 2 in equation (3.4.19), the coeffi- 
cient of a'^D'^ is positive, the coefficient of a^D^ is negative, 
and the remaining coefficient is positive. By the definition 
of A/i(a), this statement remains valid with |m| = 1 for the 
first two coefficients and for the last coefficient with an ad- 
ditional sujficient requirement a? > Cj/4. The determinant 



A of the quadratic equation (3.4.19) is positive for |m| > 2 
such that there are two positive a'^D^ (see Appendix C). 

One can further show, in a similar manner, that the 
determinant A of equation (3.4.19) is non-negative for F = 1 
and \m\ = 1 with sujficient requirements a > -\/3/2 and 
a > Ca/2. In other words, there are two positive values of 
as long as a > V3/2 and a > Ca/2. 

For F = 1, m = 0, and a a greater than a specific value 
Oc in equation (3.4.19), one can show that the coefficient of 
a*D* is positive, the coefficient of a^D^ is negative, and the 
remaining coefficient is positive. Following the same proce- 
dure of proof, one can further show that there exist two pos- 
itive roots of a^D^. For < a < Qc, however, the solution 
structure of condition (3.4.19) may have several possiblities, 
depending on the ratio = C'a/o^ . In terms of , a, and 
Mo{o), the two solutions of are 



{[l-(l-gV2)AAo(a)](a^ + l/4) 

+q^{o/' - l/4)AAo(a) + q^{a - 7/4)} 

{[l-(l-gV2)AAo(a)](aW4) 

+q\a' - l/4)AAo(a) + q^ia" - 7/4) 
-4[-2 + (a^ + l/^)NQ[a)]{q^{a^ - 1/4) 

-I 1/2T 

x[l-(l-gV2)AAo(a)]-?V4} 



x{2[-2 + (a^ + l/4)AAo(a)]}- 



(3.4.20) 



where the denominator changes from negative to positive 
between a ~ 1.79 and 1.795 approximately"*; this value of 
Q, independent of g^, is denoted by olc- We shall refer to the 
two values of as the plus- and minus-sign solutions ac- 
cording to their respective signs in front of the square root of 
the determinant in equation (3.4.20). As specific examples, 
we explored numerically three cases with decreasing values 
of ratio q^ , namely, (a) = 3.61 (Fig. la), (b) q^ — 1.0 
(Fig. lb), and (c) q^ = 0.09 (Fig. Ic). For a > Oc, there 
are indeed two positive values of as stated earlier. The 
larger and smaller ones correspond to the plus- and minus- 
sign solutions, respectively. As a — > Qc -|- 0"^, the plus-sign 
solution goes to +<x, while the minus-sign solution remains 
finite and continuous across Qc. For q^ 0, there exists a 
finite interval a; < a < q„ such that the determinant A 
is negative and there are thus no real solutions for . Ap- 
proximately, ai is between 0.595 and 0.6 while is between 
0.815 and 0.82 for = 3.61; a; is between 0.68 and 0.685 
while Q„ is between 0.92 and 0.925 for q^ = 1.0; and aj is 
between 0.96 and 0.965 while is between 1.08 and 1.085 
for q^ = 0.09. There is a systematic shift of this open interval 
{ai,au) towards larger a as g decreases. For 

OLu ^ ^ etc, 



i We used an approximate expression oijVo{a). Shu et al. (2000) 
derived a value of \ac\ = 1.759 from the same condition (a^ -|- 
l/4)Aro(ac) = 2. 
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the plus-sign solution is negative while the minus-sign so- 
lution changes from negative to positive for increasing a. 
For < a < a;, the minus-sign solution is greater than the 
plus-sign solution and the two solutions join at a = ai with 
a > 0; for = 1 and = 0.09, the smaller plus-sign so- 
lutions change from negative to positive as a increases from 
a = within this interval. All negative values of are not 
shown in Fig. 1 as they are unphysical. 

There are now two ring fragmentation regimes when a 
is sufficiently large (one for larger and one for smaller 
D^) and one modified collapse regime for sufficiently small 
a and as shown in Figure 1 for all three values. In 
the limit of q^ 0, the situation degenerates to the case 
shown in Figure 2 of Shu et al. (2000) as expected (see also 
Lemos et al. 1991 on the linear stability of axisymmetric 
scale- free disks to axisymmetric disturbances). For increas- 
ing values of q^ , the lower ring fragmentation regime asso- 
ciated with slow MHD disturbances is enlarged while the 
upper ring fragmentation regime associated with fast MHD 
disturbances is pushed upwards in the parameter scheme of 
Fig. 1. That is, one needs an even larger to access the 
upper ring fragmentation regime. 

As the MHD generalization of Toomre's Q-parameter, 
the QAf-parameter (Lou & Fan 1998a) is defined by 



we derive 



(a^ + 



2V2D{1 + I 



,2x1/2 



ttGEo 



F(l + D2_g2/2) 



(3.4.21) 



with an MSID profile (2.2.2). For axisymmetric FMDWs, 
the disk stability requires Qm > 1 (Lou & Fan 1998a). For 
the minima -D^i„ of the boundaries of the upper fragmen- 
tation regime and the relevant values of q^, the correspond- 
ing values of Qm thus obtained are fairly close to unity in 
all three cases of Figure 1. For example, with q^ — 3.61, 
D^i„ = 40.62 at a = 3.9 and Qm = 0.972; with = 1.0, 
= 15.24 at a = 3.681 and Qm = 0.992; and with 
= 0.09, L»Ln = 6.48 at a = 3.264 and Qm = 1.01. It 
appears that the critical Qm value tends to increase slightly 
with decreasing values of q^ and with decreasing values of 
critical a. This trend is also complemented by the limit- 
ing case (Shu et al. 2000) with g^ = 0, = 5.41 at 
a = 3.056, and Qm — 1.026. Thus the approximate insta- 
bility criterion Qm <1 (Fan & Lou 1996; Lou & Fan 1998a) 
appears pertinent to the upper branch of ring fragmentation 
for axisymmetric fast MHD disturbances (see Fig. 1). 

Comparing with Fig. 2 of Shu et al. (2000), the coplanar 
magnetic field modifies the collapse regime as shown in our 
Fig. 1 and introduces the lower regime of ring fragmenta- 
tion for slow MHD disturbances when a is sufficiently large. 
These instabilities do not occur when is large: > 9.2 
for = 3.61; > 3.5 for q^ = 1; > 1.2 for = 0.09; 
and > 0.932 for g^ = 0. The last result of g^ = is taken 
from subsection 4.1 of Shu et al. (2000). 

From equation (3.4.13) and 
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which is important to determine the spatial phase relation- 
ship between the surface mass density and azimuthal mag- 
netic field perturbations. In the limit of a — + 0, one has 
Z —> 0, but ill given by (3.4.22) remains nonzero. For a 
small a 7^ such that may be dropped relative to ia 
term, Z and S are phase shifted by ~ ±7r/2. 
In the limit of a — > cxd, we obtain 
BeS 



Z = 



[a - 2nGMn{a)T,or - fl^ 



(3.4.24) 



which, given stationary condition (3.4.16), can be cast into 
a convenient form to examine the spatial phase relationship 
between Z and S in an MSID plane at z = 0. 



4 PHASE RELATIONSHIPS BETWEEN 

MAGNETIC FIELD AND MASS DENSITY 

We examine spatial phase relationships among velocity dis- 
turbances, azimuthal magnetic field, and surface mass den- 
sity enhancements because they provide useful observational 
diagnostics for magnetized spiral galaxies (Beck & Hoernes 
1996; Fan & Lou 1996; Lou & Fan 1998a, 2002; Frick et 
al. 2000) . Regions of high-density gas are vulnerable to ac- 
tive star formation, while nonthermal radio-continuum emis- 
sions from gyrating relativistic cosmic-ray electrons trapped 
in a spiral galaxy would reveal regions of stronger magnetic 
field. Thus, large-scale spiral structures of optical and radio- 
continuum emissions contain valuable information of the un- 
derlying MHD (Lou & Fan 2000a, b). The mathematical 
development is somewhat lengthy and the relevant formu- 
lae are summarized in Appendices D and E for aligned and 
unaligned perturbations respectively. A reader may want to 
mainly concentrate on the flow of logics with convenient ref- 
erences to equations in the two Appendices. 

4.1 Aligned MSID Configurations 

For a full MSID of F = 1, as given by solution (3.2.10) 
depends only on two dimensionless parameters q^ = C^/a^ 
and |m|. By solutions (3.2.10), one derives equation (Dl) 
that shows a^D^ — Cj/2 > for the upper plus-sign solu- 
tion in equation (Dl). For the lower minus-sign solution in 
equation (Dl), whether inequality > C^|m|/[2(jmj — 1)] 
holds or not in the determinant A would determine whether 
a^D^ - C\/2 > or not. 

By equation (Dl), we then examine the sign of 
27rGEor = F{a'^ + a^D'^ - Cl/2) for a full MSID of F = 1. 
An addition of to the first term on the RHS of equation 
(Dl) gives equation {D2). Examine then the determinant A 
in solutions (3.2.10) or (Dl) in the form of equation (-D3). 
The sum of the last two terms in expression (-D3) is given 
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by expression (-D4). When > C4/4, we have the back- 
ground surface mass density Eo > for both the plus- and 
minus-sign solutions (3.2.10) of a^D^. For < one 
has So > for the plus-sign solution only; the minus-sign 
solution gives Eo < which is unphysical. 

To examine the spatial phase relationship between bg 
and El in the MSID plane, we use relations (3.1.8) and 
(3.1.9) derived in subsection 3.1. Specifically, we determine 
the signs of both numerator and denominator on the RHS 
of equation (3.1.9). For the numerator 2Q?r^ — C\/2 = 
2{a'^D'^ - C\/2 + Ci/4), we add Ci/4 to the first term 
on the RHS of equation (-D1) to derive equation (-D5). The 
determinant A may be cast into the form of equation (-D6). 
For - C\/4, > in equation (D6), one has both Eo > 
and 20^ D'^ - C\/2 > Q for either plus- or minus-solutions 
in equation (3.2.10). For a? — C^/4 < 0, only the plus-sign 
solution with aF'D^ — C\/A > is physically valid. 

For the sign of the denominator on the RHS of equa- 
tion (3.1.9) as given by equation {D7), a subtraction of 
Ca ^ C\/{2m^) from the first term on the RHS of so- 
lution (3.2.10) leads to expression (-D8). We now exam- 
ine the determinant A in the form of expression (-D9). 
The last term containing (3|m|'^ — 2) on the RHS of equa- 
tion (-D9) is positive for \m\ 7^ 0. Hence, the denominator 
m^ti^r^ — C\{m? — 1/2) is positive and negative for the plus- 
and minus-sign solutions in equation (3.2.10), respectively. 

To examine the sign of the background surface mass 
density Eo of a partial MSID for |m| > 1, we consider a^D^-\- 
a? — C\/2. By adding — C\/2 to the term involving B — 
A-C on the RHS of solution (3.2.12) for a^D^, one has 
inequality (DIO). In the determinant A of solution (3.2.12), 
we then consider the part as given by equation (Dll). When 

— C\/A > 0, Eo is positive for solutions (3.2.12) of both 
signs; for — C^/4 < 0, Eo is positive for the plus-sign 
solution, while Eq < for the minus-sign solution. These 
conclusions for a partial MSID (0 < F < 1) turn out exactly 
the same as those for a full MSID {F — 1). 

To examine the phase relationship between Z and S, 
equations (3.1.7) — (3.1.9) remain valid with Eo profile (2.2.2) 
for a partial MSID. Again, we examine the signs of the nu- 
merator and denominator on the RHS of (3.1.9). For the 
numerator, we consider a?D^ — C\/A by subtracting C\/A 
from the term involving B — A — C on the RHS of solution 
(3.2.12) for a^L>^. The resuhing expression is then given by 
equation (_D12). In the determinant A of solution (3.2.12), 
we then consider the part given by equation (_D13). There- 
fore, Eq and a^D^ — C\/A: are positive for both plus- and 
minus-sign solutions of (3.2.12) when — C^/4 > 0. For 

— C'^/4 < 0, only the plus-sign solution is valid with 
a^D^ - Ci/4 > 0; the minus-sign gives a < 0. 

We now examine the sign of the denominator on the 
RHS of equation (3.1.9). For this purpose, we consider 
a?D^ — [1 — 1/ {2m?)\C\. In the determinant A of solution 
(3.2.12), we consider the portion given by equation (D14), 
which vanishes for |m| = 1 and is positive for |m| > 2. 
Thus, m^f^^r^ — C\{rn^ — 1/2) is positive and negative for 
the plus- and minus-sign solutions of (3.2.12), respectively. 



Using these results in equations (3.1.8) and (3.1.9), we have 
Z and S being out of phase for the plus-sign solution of 
(3.2.12), and Z and S being in phase for the minus-sign so- 
lution of (3.2.12) when > Cl/4. When < C\/4, the 
minus-sign solution of (3.2.12) leads to a negative Eo. 

Based on these analyses and in reference to equations 
(3.1.8) and (3.1.9) for the phase relationships between Z and 
S, perturbation enhancements of be and Ei anticorrelate 
and correlate with each other for the plus- and minus-sign 
solutions, respectively, when > C^/4. When < C^/4, 
the minus-sign solution is invalid because the surface mass 
density Eo would be negative. These two distinct stationary 
aligned MSID configurations appear as results of differential 
rotation, self-gravity, and curved magnetic field. According 
to equation (3.1.4), one may write 

. _ mflrS (a^ - 27rGEor/|m| - OV^) 

Eo {Cl/2-Q^r^) ■ ^ > 

Therefore, ill oc ±rS, iR oc ±rS, and J oc —S for the 
plus- and minus-sign solutions. With these spatial phase 
relationships among stationary perturbation variables, one 
may conceive mental pictures for the two distinct types of 
stationary MSID bar configurations for |m| — 2. For the 
plus-sign case, enhancements of be and Ei are out of phase 
with each other. High-density gas regions are active in star 
formation (i.e., bright in optical and infrared bands) with 
enhanced small-scale random magnetic fields (i.e., bright in 
total nonthermal radio-continuum emissions), whereas rela- 
tively strong regular magnetic field regions should be bright 
in polarized radio-continuum emissions with higher degrees 
of polarization (less disturbed by activities associated with 
cloud and star formation on small scales). For the minus- 
sign case, enhancements of be and Ei are in phase. By the 
same rationale, the structural manifestations of bars or lop- 
sided disks in total and polarized radio emissions will be 
in competition because small-scale random motions tend to 
enhance total radio-continuum emissions while weaken po- 
larized radio-continuum emissions. Regions of strong total 
radio-continuum, optical, infrared emissions should more or 
less overlap. Depending on the level of activities in star- 
forming regions, polarized radio-continuum emissions may 
also be sufficiently strong in optically bright regions. Polar- 
ized radio-continuum emissions from relatively weak mag- 
netic field regions should show higher degrees of polarization 
because of reduced level star formation activities there. 



4.2 Unaligned Logarithmic MSID Spirals 

We first consider the important factor that appears in cri- 
terion (3.4.16) for logarithmic spirals in an MSID, as given 
by equation (-El). By equations (3.4.16) and (3.4.23), the 
sign of this factor (-E1) determines the spatial phase rela- 
tionship between be and Ei. Using the stationary criterion 
in the form of equation (3.4.19), we derive two values of 
a^D^ as given by equation {E2). The determinant A under 
the square root of solution {E2) can be shown to be non- 
negative for \m\ > 2 (see Appendix C). As fi^r^ = a^D^ , we 
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subtract Ci{l-[2(m2+a^+l/4)]"^} from solution {E2) and 
rearrange the determinant A of solution (-E2) accordingly. 
In the expression of determinant A, we consider the relevant 
part as given by equation (-E3), which is positive for |m| > 2. 
For a full MSID oi F = 1 with |m| = 1 and |mj = 0, the 
reader is referred to numerical results and discussions that 
follow equation (3.4.19). It then follows with |m| > 2 that 
the key factor (El) (see solution criterion [3.4.16]), namely 
(m^ + + l/4)r2^r^ - (m^ + - l/4)Ci , is positive and 
negative for plus- and minus-sign solutions of a?D^ in equa- 
tion {E2), respectively. 

Using dispersion relation (3.4.16) for stationary loga- 
rithmic spirals in an MSID and relation (3.4.23) between Z 
and S derived in subsection 3.4, we finally arrive at 

Z oc ±S{a + ia[C\/{2n^r'^) - 3/2]} (4.2.1) 

for the plus- and minus-sign solutions of a^D^ in equation 
{E2), respectively. For a sufficiently large a in the tight- 
winding or WKBJ regime, one may ignore the imaginary 
part in relation (4.2.1). In this regime, Z is approximately 
in-phase and out-of-phase with S for the plus- and minus- 
sign solutions, respectively. This appears to be consistent 
with the results of FMDW and SMDW analyses in the tight- 
winding approximation (Fan & Lou 1996; Lou & Fan 1998a, 
2002). For a sufficiently small a that corresponds to rela- 
tively open spiral structures, we may drop the term in 
comparison with the term proportional to ia. In this regime 
of small a 0, Z is either ahead of or lag behind S by a 
phase difference of ~ 7r/2 for both stationary fast and slow 
logarithmic spiral configurations in an MSID; this is a new 
result, not known before. 



5 DISCUSSION AND APPLICATIONS 

5.1 A discussion on bars and barred spirals 

The modal formulation and computations of Berlin et al. 
(1989a, b) and Berlin & Lin (1996) set a theoretical frame- 
work to classify morphologies of spiral galaxies. Meanwhile, 
the analytical and numerical results of Shu et al. (2000; also 
Galli et al. 2001) provide an important perspective for the 
onsets of bar-type and barred-spiral instabilities in isope- 
dically magnetized SIDs associated with nonaxisymmetric 
aligned disturbances and unaligned logarithmic spiral per- 
turbations that appear stationary in an inertial frame of 
reference. One would like to know the overall connection 
between the results and interpretations of Shu et al. (2000) 
and those of Berlin et al. (1989a, b; Berlin & Lin 1996), 
especially regarding the cubic dispersion relation and ac- 
curately solved numerical solutions in the modal scenario 
for bars and barred spiral galaxies. Among several appar- 
ent differences in the model formulations such as SIDs ver- 
sus prescribed disks, logarithmic spiral perturbations versus 
normal modes, stationarity versus time variations, exact an- 
alytical solutions versus extended WKBJ approximation as 
well as accurate numerical solutions, we single out one key 
element that, we believe, distinguishes the two independent 



lines of approach in a significant manner. The issue involves 
analyses that are somewhat technical and subtle. We shall 
describe them below. 

For SIDs and their kin (Mestel 1963), there have been 
theoretical studies on their structures and instability proper- 
ties for decades (Zang 1976; Toomre 1977; Lemos et al. 1991; 
Lynden-Bell & Lemos 1993; Syer & Tremaine 1996; Evans & 
Read 1998; Goodman & Evans 1999; Shu et al. 2000; Galh et 
al. 2001) in view of their potential applications to the struc- 
ture of lopsided or normal and barred spiral galaxies, to the 
light cusps seen in the nuclei of galaxies, and to the forma- 
tion and collapse of cloud cores in the birth of stars and plan- 
etary systems. Syer & Tremaine (1996) found semianalytic 
and numerical solutions of nonaxisymmetric stationary equi- 
libria for completely fiattened (razor-thin) power-law disks. 
The basic problem along with pertinent issues have been 
well summarized by Shu et al. (2000), and for the special 
case of index /3 = (in the notation of Syre & Tremaine 
1996), Shu et al. (2000) derived analytic solutions and cri- 
teria for both aligned and unaligned logarithmic stationary 
perturbations in razor-thin isopedically magnetized SIDs. 

For aligned perturbations, the zero-frequency solutions 
correspond to the onset of bifurcations from axisymmet- 
ric SIDs to nonaxisymmetric SIDs in close analogy to bi- 
furcations from incompressible uniformly rotating Maclau- 
rin spheroids to Dedekind ellipsoids with configuration axes 
that remain fixed in space (Chandrasekhar 1969; Binney 
& Tremaine 1987). Moreover, Shu et al. (2000) relate the 
aligned case of m = 2 to the onset of the secular barlike 
instability in the context of galactic dynamics (Hohl 1971; 
Miller et al. 1970; Kalnajs 1972; Ostriker & Peebles 1973; 
Bardeen 1975; Aoki et al. 1979; Vandervoort 1982, 1983). 

In analyses of Feldman & Lin (1973), Lau & Berlin 
(1978), Lin & Lau (1979) and Berlin et al. (1989a, b) on 
galactic density waves, there is a coefficient, usually referred 
to as B, in the standard integro-differential equations for 
spiral density waves, namely 

nr{l — i^^) r^Ki' dlnr 

(see eq. [25b] of Lin & Lau 1979 or eq. [2.5] of Berlin et al. 
1989b), where i/'^ = (lu — mQ)^ /k^. The third term on the 
RHS of equation (5.1) is related to the corotation resonance 
and the second term on the RHS of (5.1) contains the in- 
formation of the Ti or term proportional to dlnQ/dlnr 
with J2 = {2-KGfio/K^fTx (Lau & Berlin 1978; Bertin et al. 
1989b; cf. comments of Hunter 1983, Lou & Fan 1998b, and 
MYE 1999). The second term on the RHS of (5.1) may be 
split further into the form of 

^mQ.{dv/dr) _ 4m'^Q'^ dlnQ ^ Am^lu dn 
K,r{l — u^) /t^r^(l — i^^) dlnr K,'^r(l — u^) dr ' 

(5.2) 

where the first item on the RHS is the —Ti/{1 — u^) term 
(the notation Ti was introduced by Lau & Bertin [1978]), or 
equivalently, —J^/{1 — u^) term (the notation was later 
used by Lin & Lau [1979] and Bertin et al. [1989a, bj; see 
footnote 2 of Lou & Fan 1998b). To derive the cubic disper- 
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sion relation of density waves, this Ti or term plays the 
central role (Bertin et al, 1989b; Bertin & Lin 1996). The 
second term of the RHS of equation (5.2) has been some- 
how ignored and relegated to a residual R term (see equation 
[B18] in Appendix B of Lau & Bertin 1978). Without this 
parameter, the cubic dispersion relation simply reduces 
to the quadratic one in the radial wavenumber k for spiral 
density waves. In short, it is this Ti or term that is re- 
sponsible for the third small k root of the so-called "open 
mode" besides the familiar long- and short-waves (see eq. 
[3.1] of Bertin et al. 1989b). In our discussion here, this 
mechanism is referred to as the TSF effect by its physical 
nature (LYF 2001). 

The utmost reason we pinpoint this B coefficient is as 
follows. Bertin et al. (1989a, b) and Bertin & Lin (1996) 
proposed the cubic dispersion relation that was derived by 
keeping some higher-order terms in the expansion of the 
Poisson equation. In proper parameter regimes, two of the 
k roots correspond to the familiar short- and long-branches 
of spiral density waves while in others, there exists a third 
small k solution that was referred to as the "open mode" 
and was suggested to correspond to barred spiral galaxies 
on the basis of their extensive numerical calculations for 
the integro-differential equations of density waves (Bertin 
et al. 1989a, b; Bertin & Lin 1996). In contrast, Shu et al. 
(2000) studied the two cases of aligned and unaligned non- 
axisymmetric stationary perturbations in isopedically mag- 
netized razor-thin SIDs and derived the marginal criteria 
for both cases. Shu et al. (2000) suggested that the class 
of aligned stationary perturbations is related to the secular 
barlike instability in galactic dynamics. This type of bar- 
like (m = 2) instability may be suppressed by introducing 
a sufficiently massive dark-matter halo (Ostriker & Peebles 
1973). By numerically studying time-dependent overreffec- 
tion process of spiral waves across the corotation radius, Shu 
et al. (2000) further recommended that the marginal crite- 
rion for unaligned logarithmic stationary perturbations as 
an indicator for the onset of spiral instabilities in the sense 
that they be eventually amplified by a swing process near 
corotation (Goldreich & Lynden-Bell 1965; Mark 1976; Fan 
& Lou 1997). By comparing their marginal criterion (38) 
with their standard WKBJ dispersion relation (39) for spi- 
ral density waves (Lin & Shu 1966, 1968) and by identifying 
an effective wavenumber (m^ + + l/4)^/^r~^ as defined 
by their equation (40), an analogy of quadratic {rather than 
cubic) form in the effective wavenumber is apparent. Shu et 
al. (2000) proposed to use their marginal criterion (38) in 
the limit of a — > (i.e., "breathing mode" limit; Lemos et 
al. 1991) for the onset of dynamical barred-spiral instabil- 
ity. The question now is the connection or relation between 
the two seemingly different proposals for the same physi- 
cal problem in the galactic context, namely the nature of 
galactic bars and barred spiral galaxies. 

The key is to realize the following basic fact. For the 
very SID model (9) as prescribed by Shu et al. (2000), the 
third corotation resonance term on the RHS of B expression 



(5.1) vanishes.'* For their SID model (9), the second term on 
the RHS of equation (5.1) becomes 

AmQ,{dv/dr) AmQ,w . , 

«;r(l-i/2) " ~«:V2(l-i/2) ^ > 

which vanishes for uj = 0. Remarkably, a recombination of 
the two terms on the RHS of equation (5.2) makes the — Ti 
term or term needed for the cubic dispersion relation 
(Bertin et al. 1989b) disappear. In reference to equation (5.1) 
and the stationarity requirement (i.e., a; = 0), it turns out 
that B = —m?/r^. Physically, this means that the TSF is 
absent for stationary logarithmic spiral perturbations in the 
SID model (9) as given by Shu et al. (2000) and should be 
the reason that in the analysis of Shu et al. (2000), there 
is no obvious clue or counterpart for the "open mode" as a 
root of a cubic dispersion relation. In fact, equation (38) of 
Shu et al. (2000) together with their approximate expression 
(36) for the Kalnajs function give a striking quadratic form 
in terms of their effective wavenumber. 

It seems plausible that the "open mode" might corre- 
spond to the "aligned" SID case (Shu 2000, private com- 
munications). It is then the absence of the TSF for sta- 
tionary logarithmic spiral perturbations in SIDs that decou- 
ples the "aligned" and "unaligned" cases. In other words, 
the nonzero TSF for time-dependent normal mode pertur- 
bations in differentially rotating disks other than SIDs might 
somehow couple or mingle the aligned and unaligned cases 
such that some kind of dispersion relation in the spirit of 
a cubic dispersion relation emerges. Although radial propa- 
gation is absent for aligned nonaxisymmetric perturbations 
in a SID, there is an azimuthal wave propagation relative to 
the axisymmetric background SID as made clear in our anal- 
ysis. Once present, the TSF might play a nontrivial coupling 
role for nonaxisymmetric MSID perturbations. A carefully 
designed numerical test may be needed to settle this issue. 

We also note another line of reasoning. Even though 
equation (38) of Shu et al. (2000) is approximately quadratic 
in terms of their effective wavenumber k — a^Oj~^ with 
\ae\ = {m? + + 1/4)^'''^, the recommended marginal cri- 
terion (37) (equivalently, eqs. [38] or [41]) for unaligned sta- 
tionary logarithmic spiral perturbations yields only one so- 
lution of a for a given and m 7^ in a full SID (see Figure 
3 of Shu et al. 2000). In other words, by setting equal to 
a (sufficiently large, > 0.5368) constant in Figure 3 of 
Shu et al. (2000), there is only one intersection for a given 
m 7^ and the corresponding root of a for a stationary 
logarithmic spiral may be either large or small depending 
upon whether is large or small. The point here is that 
the situation of only one root of ct for stationary logarithmic 
spirals with m 7^ in a full SID does not necessarily contra- 
dict the familiar result of short and long spiral density waves 
(the two roots for the radial wavenumber k) for propagating 
WKBJ spiral patterns in disks that are not full SIDs (e.g., in 

§ We could not identify the Bs term attributed to Lin & Lau 
(1979) by Shu et al. (2000) in the paragraph following their equa- 
tion (53) (i.e. the last paragraph of their subsection 5.1). 
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partial SIDs). By the same argument, the quadratic form of 
the marginal criterion in terms of the effective wavenumber 
Qeil)"^ for stationary logarithmic spirals in SIDs might not 
necessarily contradict the cubic dispersion relation (Bertin 
et al. 1989; Bertin & Lin 1996) for nonstationary WKBJ- 
type spiral waves in disks that are not full SIDs. 

5.2 Implications for the spiral galaxy NGC 6946 

The unambiguous case of interlaced optical and magnetic 
spiral arms in the nearby spiral galaxy NGC 6946 was first 
reported by Beck & Hoernes (1996) in the almost rigidly 
rotating inner disk portion. Similar interlaced arm features 
were also suspected earlier in portions of spiral galaxies IC 
342 (Krause et al. 1989) and M83 (NGC 5236; Sukumar 
& Allen 1989). These earlier observations, the one of NGC 
6946 in particular, prompted Fan & Lou (1996) to propose 
the concepts of FMDWs and SMDWs in magnetized spiral 
galaxies (Lou & Fan 1998a; LYF 2001). Specifically, spiral 
perturbation enhancements of magnetic field and gas mass 
density of SMDWs are significantly phase shifted relative 
to each other (with a phase difference > n/2). As high- 
density gas arms are more vulnerable to cloud and star for- 
mation activities and large-scale regular magnetic fields are 
less disturbed by small-scale ISM turbulence associated with 
star formation processes along phase-shifted magnetic arms, 
this scenario naturally leads to interlaced optical and mag- 
netic spiral structures in magnetized disk galaxies. The re- 
cent wavelet analysis on multi- wavelength data of NGC 6946 
(Frick et al. 2000, 2001) revealed extended spiral arms well 
into the outer disk portion with a fiat rotation curve (e.g., 
Tacconi & Young 1989; Sofue 1996; Ferguson et al. 1998), 
which might appear to challenge our earlier proposal that 
SMDWs be largely confined within the inner disk portion of 
almost rigid rotation. 

Given the idealizations in our MSID model, the two pos- 
sibilities of stationary logarithmic spiral fast and slow MHD 
perturbations (see eq. [4.2.2]) in an MSID of fiat rotation 
curve are conceptually important for magnetized barred or 
normal spiral galaxies in general (Lou & Fan 2002) and may 
bear direct import to the multi-wavelength data analysis on 
NGC 6946 by Frick et al. (2000). In particular, for stationary 
slow MHD logarithmic spirals in an MSID with fiat rotation 
curve, spiral enhancements of magnetic field and mass den- 
sity are interlaced with a phase difference > 7r/2 by relation 
(4.2.1). For large a in the WKBJ regime, this phase differ- 
ence approaches ~ n, while for small a 7^ in the open 
regime, this phase difference approaches ~ ■k/2. Our ana- 
lytical solutions include the effects of long-range self-gravity 
and disk differential rotation. Therefore in magnetized spiral 
galaxies, slow MSID patterns can indeed give rise to radially 
extended manifestations of interlaced optical and magnetic 
spiral structures that persist well into the outer disk portion 
with a largely fiat rotation curve (Lou & Fan (2002). 

We also note by equation (4.2.1) that for stationary fast 
logarithmic spiral structures in a full or partial MSID with 
a flat rotation curve, spiral enhancements of magnetic field 



and mass density are in phase in the tight-winding regime 
but are interlaced by a phase difference of ^ 7r/2 in the open 
regime. This is a new feature that may bear consequences 
for observations of galactic structures. 

5.3 Summary 

In reference to the work of Shu et al. (2000) on stationary 
perturbation conflgurations of isopedically magnetized SID, 
we have investigated both full and partial MSID stationary 
perturbation configurations with a magnetic field coplanar 
with the disk plane. Given the model specifications, we have 
reached the following conclusions and suggestions. 

(1) For the aligned case with > Cj/4, there are two 
possible values for the rotation parameter D corresponding 
to purely azimuthal propagations of FMDWs and SMDWs, 
respectively, with distinctly different spatial phase relation- 
ships between azimuthal magnetic field and surface mass 
density enhancements; when < C^/4, there is only one 
valid value of . Also, the case of m = can be made to 
correspond to a rescaling of the axisymmetric MSID back- 
ground. Eccentric \m\ — 1 displacements may occur for un- 
constrained values in a full MSID. In the aligned case 
of barred configurations with m = 2, there are two different 
possible types of phase relationships between magnetic field 
and mass density that should be worthwhile to search for 
observationally. 

(2) For a partial SID with a fiat rotation curve, station- 
ary eccentric \m\ — 1 displacements are not allowed. For a 
partial MSID with a fiat rotation curve, stationary eccentric 
|m| = 1 displacements can no longer occur for arbitrary D^; 
they may appear only when a^D^ — C^/2. For disk galaxies, 
it is the usual case that Q^r^ = a^D^ > due to the pres- 
ence of massive dark-matter halos. One needs other physical 
conditions to produce stationary eccentric |m| = 1 displace- 
ments in disk galaxies (e.g., power-law rotation curves with 
/3 / 0; see Fig. 5 of Syer & Tremaine 1996). 

(3) For the unaligned logarithmic spiral case, there are 
two values of for |m| > 2 corresponding to FMDWs and 
SMDWs with distinct spatial pattern relationships between 
azimuthal magnetic field and mass density enhancements. 
For m — 0, there are now two ring fragmentation regimes 
and one modified collapse regime as shown in Fig. 1. In the 
absence of a coplanar magnetic field, the small-D^ ring frag- 
mentation regime disappears. In terms of phase relationships 
between azimuthal magnetic field and surface mass density 
enhancements, the results in full and partial MSIDs are ba- 
sically the same. That is, in the tight- winding regime, the 
two enhancements are in phase for fast MSID configurations 
but are out of phase for slow MSID configurations. For open 
structures, the two enhancements are phase shifted by ~ 7r/2 
for either fast and slow MSID configurations. 

(4) The MSID virial theorem (3.3.2) is suggestive that 
the modified ratio T/\W - M\ = D'^/[2{1 + D^)] be cru- 
cial to determine the stability property of an MSID with a 
coplanar azimuthal magnetic field, where T > is the ro- 
tational kinetic energy, W < is the gravitational potential 
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energy, and > is the magnetic energy of the entire 
MSID system. 

(5) Regarding the conceptual connection between the 
cubic dispersion relation in the modal formulation (Bertin 
et al. 1989a, b; Bertin & Lin 1996) and the perspective of 
stationary SID configurations (Shu et al. 2000), we suggest 
that the absence of the tangential shear force (TSF) in a SID 
of a fiat rotation curve decouples the bar modes (i.e., aligned 
configurations) from the spiral modes (i.e., unaligned loga- 
rithmic spiral configurations). This property is also carried 
over to our investigation of stationary MSID configurations 
with a coplanar magnetic field. 

(6) While our model formulation is idealized, the results 
provide a conceptual basis and useful clues for diagnostics 
of galactic bars and lopsided or barred and normal magne- 
tized spiral galaxies. For example, for stationary unaligned 
logarithmic spiral patterns of SMDWs, the interlaced opti- 
cal and radio-continuum spiral structures may well extend 
into the disk domain with a largely flat rotation curve as in 
the case of NGC 6946 (Lou & Fan 2002). 
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6 APPENDIX A 

The proof of determinant A > in equation (3.2.10) goes as 
follows. Before factorization (3.2.8) we solve equation (3.2.7) 
with F = 1 for two values of a^D^ as 

a^D^ = {B-A-C± A'/2)/[2(2 -m' - \m\)] , (Al) 

where 

A = m^[a^ - a^/\m\ + Ci/(2]m])] > , (A2) 

B = 3Ci/2 > , (AS) 

C = Ci(m' -l/2)(l + l/|m|) >0 , (A4) 
and the determinant A 

A = {A + B-Cf + C\{rr? + \m\ - 2) 

-l-2CA(l-2/m^)^>0 (A5) 
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for |m| > 1. In solution (Al) for a?D^, one can factor out 
(|m| — 1) based on the following three factorizations, 

[1 - l/(2m^)]^ - Ci/4 = (|m| - l)[(|m| 4- l)^/m^ 

+ {a - C\/2)/{2\m\)] , {A&) 

2(2 - |m|^ - |r7i.|) = -2(|m| - l)(|m| + 2) , (A7) 

A-B+C = {\rn\ - l)[\m\a^ + C\/2 

+Ci(|m| + lf/\m\ - Cl/(2|m|)] . (A8) 

This then completes the proof that the determinant A is 
non-negative for \ni\ > 1 in solution (3.2.10) for a?D^ . 

7 APPENDIX B 

As in Appendix A, we introduce handy notations 

A = rn\a - aFl\m\ + C\Fl{2\m\)\ > , (Bl) 

B = 3Ci/2>0, {B2) 

C = C\{m^ -l/2){l + F/\m\)>Q . {B3) 
The determinant A in solution (3.2.12) may be written as 

A={A + B-Cf + C\{m' + F\rn\ - 2) 

+2Cl(l-2/m^)^>0 (B4) 

for \rn\ > 1. Both solutions of a^D^ in (3.2.12) are positive. 

8 APPENDIX C 

Parallel to the proofs given in Appendices A and B, we out- 
line the proof here by introducing a set of handy notations 

A=[a^- (a - Cl/2)FAfmia)]{m^ + a' + 1/4) , (CI) 

B = C\{n? +c? - l/4)FAA^(a) , (C2) 

C = Cl(m^-f Q^-7/4) . (C3) 
The determinant A of equation (3.4.19) is given by 

A = (^ + S + C)^ - 4[m2 - 2 + (mV l/4)FAA„(a)] 
X \C\{m^ + - \l\)AI(m + a V 1/4) - Cl/4] > (C4) 

for \m\ > 2. One needs to shuffle and regroup a few terms by 
noting a -4AB term and adding 2C{A-B) and -2C{A-B) 
terms and so forth. It is useful to refer to Figure 1 of Shu et 
al. (2000) for the magnitudes of J^m{a) (m = 0, 1, 2). 

9 APPENDIX D 

For aligned perturbation conflgurations constructed from a 
fuU MSID of F = 1, as given by equation (3.2.10) de- 
pends on = Ca/cl^ and \m\. By solutions (3.2.10), 

a^D^ - Cl/2 = [(|m| + l)C'i/2 + a^\ni\'^]/[2\m\{2 + \m\)] 

±[2|7Ti|(2 + \Tn\)]-^{[Cl{\m\ + l)/2 + 

-2C'\\m\{2 + |m|)(|m| + l)[a^(|m| - 1) - Cl\m\/2]}^^'^ 
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An addition of to the first term on the RHS of solution 
(-D1) gives 

[(|m| + l)Ci/2 + a^|m|^]/[2|m|(2 + |m|)] + 

= [(|m| + l)Ci/2 + 3a^|m|^+4a^|m|]/[2|m|(2 + lml)] . (D2) 

Examine then the determinant A in solutions (3.2.10) or 
(Dl) in the form of 

[(jm| + l)Cl/2 + 3a^\mf + 4a^\m\f 

-^a(rn + 2|m|)(|m| + X){2a\'m\ + C\/2) 

-2C\\m\{2 + |m|)(|m| + l)[a^(|m| - 1) - C\\m\/2] . (D3) 

The sum of the last two terms of expression (-D3) is 

-d,{\Tri\^ +2\ni\){\m\ + l)\'m\{a -C\/4){a +C\/2) . (D4) 

When a? > the background surface mass density 

Eo > for both the plus- and minus-sign solutions (3.2.10) 
of a^D^ For < Ci/4, one has Eo > only for the plus- 
sign solution, while the minus-sign solution gives Eo < 0. 

For the numerator 2!^V^ - Ci/2 = 2(a^D^ - Cl/2 + 
Ci/4) on the RHS of equation (3.1.9), we add Ci/4 to the 
first term on the RHS of solution (Dl) to derive 

(|m|-hl)Ci/2-ha'|m|2 Cl 



+ 



2|m](2-^|m|) ' 4 

_ (Imp + 3imj + l)Cl/2 + a^\rn\^ 



> 



(D5) 



2|mj(2 + \m\) 

The determinant A may be rearranged into the form of 

A = [(Imp + 3|m| + l)Ci/2 + a'^\m\^f 

-Ci(|mp + 2|m|)(3|mp - 2){a^ - Ci/4) . (D6) 

For — Ca/4 > in equation (-D6), one has both Eo > 
and 2o^D^ — C\/2 > for either plus- or minus-solutions 
in equation (3.2.10). For a? — (7^/4 < 0, only the plus-sign 
solution with a?D^ — C^/4 > is physically valid. 

The denominator on the RHS of equation (3.1.9) is 

m^J^V^ - C\(m^ - 1/2) = m^[nV^ -C\ + C\/[2m^)] . 

iD7) 

A subtraction of C\ — C\/ {2m?) from the first term on the 
RHS of solution (3.2.10) gives 

C\/2 + 5Cijm|/2 + (a" + Ci)|mp 



2|m|(2 + |m|) 
3Ci(l - |m|)/2 + {a' - Ci)|mp -I- 2Ci/lm 



(D8) 



2|ml(2 + |mj) 
The determinant A is now examined in the form of 

A = [3Ci(l - |m|)/2 + {a" - Ci)|mp + 2C\/\m\f 

+C\{\fn\ + 2){\m\ + l)(3|mp - 2)/|mp > (D9) 

for m 7^ 0. 

For the sign of the background surface mass density Eq 
of a partial MSID with |m| > 1, we consider oF'D^ + — 
C\/2. By adding a? — C\/2 to the term containing B — A — C 
on the RHS of solution (3.2.12) for a^D^ , one has 

m^a^ + 2a (jr? -f F|m|/2 - 2) (|mj - \m\-^)C\F 12 



2(m2 -hF|m| - 2) 



> . 
(DIO) 



In the determinant A of solution (3.2.12), we then consider 
the following part 

-2{a - C\l2)(Tn + F\-m\-2) (DU) 
X [2(a^ - Ci/2)(m' + F|m| - 2) + 2(^ - S + C)] 
-f4Ci[(m^ - \l2)Alm^ - Ci/4](2 - m^ - F\m\) 



= -8 m 



+ F\m\ - 2)(m^ - \){a + C\l2){a - Ci/4) 



where A, B, C are defined by equations (-B1), {B2), (-B3) 
in Appendix B. When — C\/A > 0, Eo is positive for 
solutions (3.2.12) of both signs; for a? — C^/4 < 0, Eo is 
positive for the plus-sign solution, while Eo < for the 
minus-sign solution. 

For the numerator on the RHS of (3.1.9), we consider 
o?D^ — C\/A by subtracting C^/4 from the term involving 
S - ^ - C on the RHS of solution (3.2.12) and obtain 
a^im^ - \ni\F) + (m^ + 2\m\F - 2 - F/\m\)C\/2 



2{m? + F\m\ - 2) 



> . 
{D12) 

For the determinant A of solution (3.2.12), we then consider 
the following part 

C\{m^ + F\m\ ~ 2)[A ~B+C- Cl(m^ + F\m\ - 2)/4] 
+AC\[[rn - l/2)A/rn - Ca/4](2 - rn - F\m\) (D13) 
= -~C\{rn + F\m\ - 2)(3m^ - 2)(a - C\/'i){l - F/\m\), 

where A, B, C are defined in Appendix B. 

For the sign of the denominator on the RHS of equa- 
tion (3.1.9), we consider a'^D^ - [1 - l/(2m^)]Ci. In the 
determinant A of solution (3.2.12), we consider the portion 

(2 - l/m^)Ci(m^ + F\rn\ - 2) 
X [2{A -B + C)-{2- \lm^)C\(m + F\Tn\ - 2)\ 
+4Ci[(m^ - \l2)Alm^ - Ci/4](2 - m^ - F\m\) 



= C\(jn ^F\m\ - 2)(jn - l)(3m^ - 2)/r: 



(Z?14) 



Thus, rr^^r^ — C\(rn^ — 1/2) is positive and negative for 
the plus- and minus-sign solutions of (3.2.12), respectively. 



10 APPENDIX E 

We consider the following factor that appears in criterion 
(3.4.16) for logarithmic spirals in an MSID, 



{rn + a + 1/4)^ - (m^ + c 



l/4)Cl 
1 



= (m^ + + l/4)<^ n r - 1 — - , 

^ ' ^\ [ 2{vn? + a^ ^ 1/4) 

Using the stationary criterion in the form of equation 
(3.4.19), we derive two values of a^D^ 

2p2_ ^ + B + C±A^/^ 

" 2[m2-2 + (m2 + Q2 + l/4)FAA™(Q)] ' ^ ' 

where A, B,C, and A are defined by equations (CI) — (C4) in 
Appendix C. The determinant A is non-negative for \m\ > 2. 
We subtract Ci{l - [2(m^ + a'^ + 1/4)]"^} from solution 
{E2) and rearrange the determinant A of solution {E2). In 
the determinant A, we consider the following part, 

[m^ -2+ (m' + + l/4)FAA™(a)] 
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Figure 1. Parameter regimes in terms of D^, a, and = C'^/a^ 
separated by stationary MSID configurations for unaligned log- 
arithmic spirals with m = 0. (a) = 3.61; (b) q^ = 1.0; (c) 
g2 = 0.09. 

x[2-l/(m2 + a2 + l/4)]Ci 

x{2{A + e + C) - [m'' - 2 + (m' + l/4)FAA™(a)] 
x[2 - l/(m^ + a' + l/4)]Cl} 
-4[m^ - 2 + (m" + + l/4)FMn(a)] 
x[Ci(m' + a' - l/4)yi/(m" + a" + 1/4) - Cl/4] 
= Clfm" - 2 + (m^ + a" + l/4)FAA„(a)] 
x{[2 - l/(m^ + 1/4)] [2(a^ + 1/4) 

+ (m^ - 2)/(m2 + + 1/4)] + 1} , (F3) 

which is positive for |m| > 2. 
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